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I. INTRODUCTION 

The study of the emergence of classical features in systems 
ruled by quantum mechanics is as old as quantum mechan- 
ics itself. When the quantum system is isolated and the evo- 
lution unitary, these features appear in the WKB semiclassi- 
cal limit, which is of paramount importance in establishing 
the quantum classical correspondence in integrable systems. 
In its more modern form, the EBK quantization rule [j]], it 
shows the direct connection between tori in phase space and 
quantized eigenfunctions in Hilbert space. In chaotic systems, 
the relationship is more subtle and is embodied in the cele- 
brated Gutzwiller trace formula |2], relating sets of unstable 
periodic orbits to the density of states. The limits of appli- 
cability of these semiclassical methods and the insight they 
provide on the quantum dynamics of isolated chaotic systems 
has inspired most of the recent research in the area of quantum 
chaos. 

In open quantum systems, on the contrary, the emergence 
of classical features has been studied mainly in the time evo- 
lution of a different set of observables, most notably the rate 
of linear entropy growth (or purity decay) |3j, the Loschmidt 
echo or fidelity El 01 and the decay of correlations Q] . 
These studies have demonstrated that in certain well defined 
regimes for chaotic systems the classical Lyapunov exponent 
governs these rates and that the evolution of localized quan- 
tum densities in phase space becomes classical. 

In this article we consider this question from the point of 
view of the spectral properties of the classical and quantum 
propagators. Classical densities evolve according to the Li- 
ouville equation whose solution can be written in terms of 
a propagator called Perron-Frobenius operator(PF)|41J. It is 
unitary on L 2 . However, for chaotic systems, correlation func- 
tions exhibit oscillations and exponential decay. The decay 
rates are given by the poles of the resolvent of the PF opera- 
tor, the so-called Ruelle-Pollicott (RP) resonances|9]. By lim- 
iting the resolution of the functional space, one can effectively 
truncate the PF to a nonunitary operator of finite size (say 
TV x N) with a spectrum lying entirely inside the unit circle, 
except for the simply degenerate eigenvalue 1 . In the, prop- 



erly taken, limit of infinite size and no coarse graining, the 
isolated eigenvalues turn out to be the RP resonances 1 171 Hal . 
As shown in |8], the linear entropy and the Loschmidt echo, 
for asymptotic times much longer than the Ehrenfest time 1 , 
also show characteristic decay rates governed by the classical 
RP resonances|9|. Experimental evidence of this dependence 
on RP resonances was observed for the first time in 1 10]. Our 
approach is similar in spirit to the calculations performed on 
the sphere for the dissipative kicked top by Haake and col- 
laborators flllfl2lfl3ll . Fishman fl4lll5ll and, for the baker's 
map, by Hasegawa and Saphir 1 16]. We model the unitary dy- 
namics by means of a quantum map and implement a diffusive 
superoperator represented by a Kraus sum. Two recent works 
by Blank, et al. | 17| and Nonnenmacher| 18| provide a rigor- 
ous theoretical underpinning to our calculations for quantum 
and classical maps on the torus. 



The plan of the paper is as follows. Section II provides a 
short account of the quantization procedure for maps acting on 
a classical surface with periodic boundary conditions in both 
coordinates and momenta, i.e. a torus. In Sec.III we imple- 
ment the open system dynamics with the definition of a diffu- 
sion superoperator represented as a Kraus sum. The general 
spectral properties of both the unitary and the noisy part, as 
well as those of the combined action are studied. Sec. IV deals 
with the relationship between the classical and the quantum 
resonances and, utilizing recently proved theorems il7l fl8il . 
how they coincide in specific ranges of h and of the noise 
strength. As a consequence we show that the asymptotic time 
behaviour of several quantities is classical and depends on the 
Ruelle-Pollicott resonaces closer to the unit circle. A numeri- 
cal method that allows the calculation of the leading spectrum 
of resonances is developed. Sec.V illustrates this correspon- 
dence taking the perturbed Arnold cat map as an example. We 
relegate to the Appendix some notation concerning the spec- 
tral decomposition of superoperators and the details of the nu- 
merical method. 
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1 In our case the Ehrenfest time rig is related to the time it takes for an 
initially localazed package to reach the borders of phase space due to ex- 
ponential instability (it is sometimes called "log time"). 
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II. UNITARY DYNAMICS ON THE TORUS T 

We picture the classical phase space as a square of unit 
area with sides identified. The classical transformations will 
map this square onto itself, thus providing a simple model 
of Hamiltonian area preserving dynamics. The fact that the 
phase space has finite area brings some well known special 
features to the quantization that we briefly review. Refs. 
lll9l.l2"oTl provide a more extensive account. 



A. The Hilbert space 

As the phase space has finite area, which we normalize to 
unity, the Hilbert space Hn is finite and its dimension N sets 
the value of Planck's constant to h = (2ttN)^ 1 . The po- 
sition and momentum bases are then sets of discrete states 
| q), | p ), q,p = 0, . JV — 1 which are related by the discrete 
Fourier transform (DFT) of dimension N 



[p\q) 



i 



-pq 



N 



(1) 



A vector | <p(t) ) in Hn characterizes pure states of the sys- 
tem and can be represented by the amplitudes (q\<p), (p\f) in 
the coordinate or momentum basis, respectively. 

In the description of open systems it is imperative to rep- 
resent states by a density operator p. They form a subset of 
self-adjoint, positive semidefinite matrices with unit trace in 

def 

Hn 2 = Hn <8> H* N , the space of complex N x N matrices, 
usually called in this context Liouville space. While Hilbert 
space is the natural arena for unitary dynamics, this much 
larger Liouville space sets the stage for the more general de- 
scription of open quantum dynamics. It acquires the structure 
of a Hilbert space with the usual introduction of the matrix 
scalar product 



(A,B) — Tt(A^B). 



(2) 



where A,B£ Hn 2 - Linear transformations in this space are 
termed superoperators, they map operators into operators and 
are represented by N 2 x N 2 matrices. In Appendix A we re- 
view the various notations and properties related to this space. 



B. Translations on the torus 



The usual translation operator in the infinite plane K 2 is 



= e --k<iP s -kpQ e 4i 
= U q V p e^ qp , 



QP 



(3) 



(4) 



where U and V generate shifts in the position and momen- 
tum eigenbasis respectively. On the torus the main differ- 
ence is that the infinitesimal translation operators P, Q with 



the usual commutation rules cannot be defined because posi- 
tion and momentum eigenstates are discrete. However, finite 
translation operators U and V that have the property 



V P U" = U q V p e. i ^~ qp 



(5) 



can be defined and they generate finite cyclic shifts in the re- 
spective bases ll2lll . The N x N grid of coordinate and mo- 
mentum states constitutes the quantum phase space for the 
torus. Eq. Q allows a definition of a translation operator 
T(q tP ) ■ Hn — * Hn, q,p integers, analog to Eq. 10. The 
action of T on position and momentum eigenstates is 



T (qi,Pl)\l) 

T( quPl )\p) 



exp 
exp 
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P + Pi). (7) 



These equations confirm that T( ?liPl ) are indeed phase space 
translations. They satisfy the Weyl group composition rule 



^ 1 (9i,pi)^(g2,P2) 



Ti 
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(8) 



The iV 2 translations T( q p ),p,q = 0, ..JV — 1 satisfy the or- 
thogonality relation 



(9) 



thus constituting an orthogonal basis for the Liouville space 
Hn 2 - 

The expansion of any operator A in this basis constitutes 
the chord 1 19] or characteristic function representation. This 
representation assigns to every A £ Hn 2 the c-number func- 
tion a(q,p) = N~ x Tx{AT, ,) and therefore every operator 
has the expansion 



A = J2a(q,p)f (q 



p)- 



(10) 



For representation purposes we also use a basis of "phase 
point " operators that constitute the Weyl, or center II 9II. rep- 
resentation. In this basis the density operator is the discrete 
Wigner function of the quantum state. The peculiar features 
of the discrete Wigner function for Hilbe rt sp aces of finite di- 
mension have been described recently in [20]. 



C. Unitary Dynamics: quantum maps 

A classical map is a dynamical system that usually, but not 
exclusively, arises form the discretization of a continuous time 
system (by means of a Poincare section, for example). Al- 
though it is always possible, by integration of the equations 
of motion, to derive the map from a Hamiltonian, this con- 
nection is rather involved and in many instances it is more 
useful to model specific features of Hamiltonian dynamics by 
directly specifying the map equations without going through 
the integration step. The same is true in quantum mechanics: 



3 



instead of modeling the Hamiltonian operator and integrating 
it to obtain the unitary propagator, it is simpler to model di- 
rectly the unitary map. Classically an area preserving map is 
characterized by a finite canonical transformation and the cor- 
responding quantum map is the unitary propagator that rep- 
resents this canonical transformation. There are no exact and 
systematic procedures to realize this correspondence. On the 
2-dimensional plane M 2 relatively standard procedures (see 
l22lo give an approximation of the propagator in the semi- 
classical limit as 



U(q u q2)=[ r 



d 2 s 

h dqidq^ 



1/2 



cxp 



■zS(q 2 ,qi) 
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where S is the action along the unique classical path from qy 
to q-2, and where for simplicity we do not consider the ex- 
istence of multiple branches and Maslov indices. Only for 
linear symplectic maps on M. 2 this unitary propagator is ex- 
act, and then S is minus the quadratic generating function of 
the linear transformation. On the other hand, several ad-hoc 
procedures for the quantization of specific maps have been 
devised: some integrable (translations]^]! and shears) and 
chaotic maps, such as cat maps l23ll24ll . baker maps f25ll26ll 
and the standard map l27ll . Also all maps of a "kicked" na- 
ture, realized as compositions of non-commuting nonlinear 
shears can be quantized, as well as periodic time dependent 
Hamiltonians I28I1 . Once the quantum propagator has been 
constructed, the advantages of using quantum maps to model 
specific features of quantum dynamics become apparent. The 
propagator U is a unitary NxN matrix, propagation of a pure 
state is achieved simply by matrix multiplication, and finally 
the classical limit is obtained by letting N — > oo. 

In Liouville space the evolution of the density operator p by 
the map U is given by 

p = UpU f . (12) 

As a linear map acting on 7i W 2 Eq. (I12> can be written as 

p' = U®UHp) = U(p), (13) 

In what follows the notation U ® W is meant to be equivalent 
to the Ad(U) notation customary in group theory. The linear 

operator U = f U ® is a unitary iV 2 x N 2 matrix. 



m. NOISY DYNAMICS 

Realistic quantum processes always involve a certain de- 
gree of interaction between system and environment. In this 
case the evolution of the system is not unitary and requires 
a description in Liouville space. This loss of unitarity leads 
to decoherence and to the emergence of classical features 
in the evolution. When the environment is taken into 
account the evolution of the system is governed by a mas- 
ter equation, which takes the form of a hierarchy of integro- 
differential equations. A drastic simplification follows from 
the assumption that the environment reacts to the system suf- 
ficiently fast, in such a manner that the system looses all prior 



memory of its state, i.e. that the evolution is Markovian. 2 . 
The resulting Lindblad equation 1 3 3l 13411 

i - 

3 

determines the evolution of open quantum systems through a 
Hamiltonian H that governs the unitary noiseless evolution 
and the Lindblad Li operators that model the interaction with 
the environment. The particular structure of the equation en- 
sures that the evolution preserves the total probability, the pos- 
itive semi-definiteness and hermiticity of the density matrix. 
The infinitesimal propagator is a linear operator in Liouville 
space which can be integrated to yield a finite linear mapping 



S(/5o). 



(15) 



This mapping, as a reflection of the structure of the Lindblad 
equation, also has a particular form that guarantees the preser- 
vation of the general properties of the density operator. The 
general form, called the Kraus representation |35] is 



OS)- 



(16) 



The only further restriction on the M M operators arises from 
the preservation of the trace that requires 



E 



I. 



(17) 



In what follows we will select them from a certain complete 
family, with a specific norm. In that case the representation 
takes a more general form 



S(p) = J2^M»pMt 



(18) 



where now the positivity requirements are c M > and 

Y^^MlM^I- (19) 

Within this general framework, just as in the case of quantum 
maps, we have the choice of modeling the noise through the 
Lindblad operators or directly in terms of the integrated form 
via the Kraus operators. In what follows we choose the latter 
and thus we model the evolution by specifying a quantum map 
to represent the unitary evolution followed by a noisy step, 
modeled by its Kraus superoperator form. An evolution of the 
density matrix specified in this way is known in the literature 
l36U37ll38u39ll as a quantum operation. It includes the special 
case of unitary evolution when the sum is limited to only one 
term. In that case, and only then, the dynamics is reversible. 
A general superoperator has no inverse. 



2 For a detailed description of quantum noise and quantum Markov processes 
see OllESl- 
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A. Quantum coarse grained dynamics 

We are interested in modeling the effect of a small amount 
of noise on the evolution of an otherwise unitary quantum map 
l& [o[l&l4fj|| . We assume that the one step propagator results 
from the composition of two superoperators. The first is the 
unitary propagator U and the second is a quantum diffusion 
superoperator D e , defined by 



E 



Ce{q,p)f( q , P ) ®T (qip)> 



(20) 



that introduces decoherence. The linear form of the full prop 
agator is 

L = D £ oU 



and its action on a density matrix p is 

U(/5) = ^c e (q,p)f (q<p) UpU' ! f (llp y 



(21) 



(22) 



As the Kraus operators in this case are unitary the condition 
( 1191 becomes simply 



<-(q,p) = !• 



(23) 



Subject to this condition c e (q,p) can be an arbitrary positive 
function of q and p. Its significance in terms of coarse graining 
is clear: as long as c e (q,p) is peaked around (0,0) and of 
width e, the action that follows the unitary step consists in 
displacing the state incoherently over a phase space region of 
order e. To avoid a net drift in any particular direction, c e (q, p) 
must be an even function of the arguments q and p. From this 
imposition and the fact that 



T. 



'■(q,p) ~ T (-1-P) 

it follows, from the properties of the matrix scalar product 
(0 that D e is hermitian (see App.lXlfor details on the scalar 
product). 

The spectral properties of the separate superoperators U and 
D e are simple to obtain. If the Floquet spectrum of the quan- 
tum map is 

U\ 4>k ) = | fc ) (25) 
then the spectrum of U is unitary and given by 

U(|^)(0 J -|) = e l «*-^|^)(^|. (26) 

To obtain the spectrum of D e we use the composition rule 
to show that 



(24) 



^(9,P)^W)^(g,p) 



exp 



,2tt 

1 77 



(vq - up) 



(27) 



We then derive 



^C e (q,p)f (q:P) f M fl qp) 
q,p 

^c e (g,p) exp 

9,P 



T, 



(28) 



Therefore, the A^ 2 eigenvalues of D e are given by the 2D dis- 
crete Fourier transform c e (fi,u) of the coefficients c e (q,p). 
The eigenfunctions are the translation operators themselves. 
Hence, using the bra-ket notation described in appendix 1X1 
the spectral decomposition of D e is 



(29) 



in analogy with Eq. JA10> . 

Physically, the action of D c is quite simple in the chord 
representation dlOl : if p is expanded as 



then 



°c(p) = 'Yj C z{P'i v )PwT( 



(30) 



(31) 



Thus the coefficients in the chord representation are sup- 
pressed selectively according to c e (p,, v). 





FIG. 1: The left pane c t (p,,v) shows the eigenvalues of D e for 
e — 0.15 and N = 100. The DFT of this function generates the 
coefficients c e (q, p) (right pane) of the Kraus representation of D e of 
Eq. J2Q). 



It is evident from Eqs. (I20> and \29\ that the diffusion su- 
peroperator thus defined can be specified indistinctly either by 
c e (q,p) or by c e (/i, v) . For an efficient numerical implemen- 
tation of its action we have found convenient to specify the 
latter as 



+sin [ITU 



/N]) 



(32) 



This is a smooth Gaussian-like periodic function of the integer 
variables p and v. For large values of N it is very close to the 
Gaussian 



(33) 



This means that the action of D e will leave essentially un- 
altered the coefficients p^ v in a region of size ~ l/(eJV) 
(Fig. [2 left) around the origin while strongly suppressing 
those outside. The backward DFT of c e (p, v) does not have a 
simple analytic expression but from general properties of the 
DFT it will also be a Gaussian like function with the comple- 
mentary width ~ e/2-7r (Fig. fright). 

The action of D e progressively washes out the quantum in- 
terference. This fact is clearly seen if the density matrix is 
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(a) (b) 




FIG. 2: Display of the action of D e . Panel (a) shows the Wigner 
function after the step U has been applied to a position state (po = 
I QQ ){lo D- Panel (b) shows the state after the full propagator L s = 
D e ol) has acted. The map is the perturbed cat map of Eq. I64l fc = 

0.02, N = 60, e = 0.25. 



represented by the Wigner function. On the torus the Wigner 
function exhibits two different types of interference. The 
stretching and folding produce quantum interference between 
different parts of the extended state. Additionally, the peri- 
odicity of the torus introduces interference between the state 
and its images. In Fig. [2] we show the difference between the 
unitary and the noisy evolution of a coordinate state by a non- 
linear map. The two types of interference are clearly seen. 
The long wavelength fringes on the convex side are produced 
by nonlinearities. The short wavelength fringes correspond to 
the images. The effect of the noise can be seen on (Fig.[3b)) 
: the classical part of the state (in white) has been broadened 
and the long wavelength interference has been significantly 
erased. This process continues at each step of the propaga- 
tion and the quantum state becomes more and more mixed 
and more and more classical. 



In the limit of large N we can thus write 

fNe\ 2 , , e 2 
-{-) ^M<- T (35) 

The singular values accumulate near the origin, thus 
forcing most of the eigenvalues of L e to be near zero. 
On the other hand the allowed eigenvalue region ex- 
tends exponentially close to the unit circle in the limit 
0. 

• The superoperator is not normal, and therefore has dis- 
tinct left and right eigenoperators corresponding to each 
eigenvalue. The left and right eigenvalue problems are 
then posed as follows for each pair of complex conju- 
gate eigenvalues A, A* 

L e R i = \ i R i L e AJ = A| J Rj (36) 
U e L i = X* i L i L\L\ = \L\ (37) 

where Li , Ri conform a biorthogonal set 
Tr (L\R 3 ) = Tr(L^) = Tr (l\r\) = 5^. (38) 

and we assume the normalization Tr ^l\L^J = 
Tr (^R]Rij = 1 In particular, corresponding to A = 1 

we choose Lq = Rq = I/N and therefore all the re- 
maining eigenoperators are traceless. 

• The spectral decomposition of L e then becomes 

L e = '£\Ri)\i(L i \. (39) 



B. Spectrum of the quantum coarse gained propagator 

In this section we study the general features of the spec- 
trum of the combined action of the unitary map and the coarse 
graining operator, given by (Eq. i22\ ). For finite values of e 
and N, L e is a convex sum of unitary matrices and is therefore 
a completely positive, contracting superoperator. Its spec- 
trum has the following properties: 

• It is unital i.e. it has a trivial non-degenerate unit eigen- 
value corresponding to the uniform density = I/N. 

• The remaining spectrum is entirely contained inside the 
unit circle and symmetric with respect to the real axis. 
The pair of complex conjugate eigenvalues correspond 
to hermitian conjugate eigenoperators. 



• As y (L e oLj) = D e the eigenvalues of D e are also 
the singular values of L £ . Therefore the spectrum is 
contained exactly in the annulus 

e -(^) 2 <|A,|<e-H^->/"]) (34) 



The exact numerical calculation of the spectrum is ham- 
pered by the need to diagonalize very large non-hermitian ma- 
trices of dimension N 2 x N 2 for values of N large enough to 
extract semiclassical features from the spectrum. In Section 
II V Bl we develop a method, specially adapted to chaotic sys- 
tems, that takes account of the dynamics of the map to extract 
the part of the leading spectrum relevant to asymptotic time 
behavior. 



IV. QUANTUM CLASSICAL CORRESPONDENCE 

Chaotic evolution in phase space implies exponential 
stretching and squeezing of initially localized densities. On 
a timescale of the order of the Ehrenfest time tf%, significant 
quantum corrections to the classical evolution inevitably ap- 
pear. However, essentially classical features emerge from 
quantum chaotic dynamics when decoherence is introduced, 
even in the limit of no decoherence. In this section we relate 
the spectra of the propagators of densities (both classical and 
quantum) with the underlying, mainly asymptotic, behavior of 
time dependent quantities. 
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Consider the classical analog for the propagation of densi- 
ties in phase space. If /(x) is a classical map, and x = (q,p) 
a point in phase space, then the evolution of a probability den- 
sity is governed by 

P'(y) = [S(y- /(x))p(x)dx - [Cp] (y) (40) 

where y — (q',p') and C is the Perron-Frobenius (PF) 
operator| 4 1J. It is unitary on the space of square integrable 
functions L 2 , and infinite dimensional. However, one is 
mostly interested in the decay properties of observables much 
smoother than L 2 . When the functional space on which C op- 
erates is restricted by smoothness, the spectrum of PF changes 
drastically, moving to the inside of the unit circle. This 
smoothing can be attained by convolution with a self adjoint 
compact (on L 2 ) coarse graining operator T> e I18l 14211 . where e 
is the coarse graining parameter. The coarse grained PF takes 
the form 

C e = V e o C, (41) 

(notice the analogy with Eq. i2H ). T> £ damps high frequency 
modes in L 2 and thus effectively truncates £ to a nonunitary 
operator. There is substantial difference, however, between 
the spectrum of the PF for a regular map and for a chaotic 
map. As the coarse graining e tends to zero, parts of the spec- 
trum of C e for a regular map can be arbitrarily close to the 
unit circle. On the contrary, for a chaotic map there is a finite 
gap for any value of e > 0. The isolated eigenvalues which 
remain inside the unit circle as e — > are the Ruelle-Pollicot 
resonances. Rugh|43] and more recently Blank, etal. ifTvtl 
made formal descriptions of the spectrum of PF for Anosov 
maps on the torus using tailor-made Banach spaces adapted to 
the dynamics. Moreover Blank, et al. use this to analyze res- 
onances of noisy propagators and prove that these resonances 
are stable, i.e. independent of e in the limit of small coarse- 
graining . Blum and Agam|44] proposed a numerical method 
to approximate the classical spectrum using similar concepts. 

A formal and very thorough recent work by 
Nonnenmacher| 18| explores the characteristics of prop- 
agators, both classical and quantum, with noise for maps on 
the torus, both regular and chaotic. In that work its is proved 
that, in the limit N — > oo, the spectrum of the coarse grained 
quantum propagator L e , for fixed e, tends to that of the 
coarse grained PF C t {\ 18] Theorem 1). These two theorems, 
taken together, provide a solid framework for the numerical 
calculation of quantum resonances of torus maps and of their 
classical manifestations. 



A. Asymptotic behavior 

The time evolution the von Neumann entropy was used by 
Zurek and Paz|3] to characterize quantum chaotic systems. 
They conjectured that the rate of increase of the von Neu- 
mann entropy of the decohering (chaotic) system is indepen- 
dent of the strength of the coupling to the environment and is 
ruled by the Lyapunov exponents. Thus classicality emerges 



naturally and correspondence even for chaotic systems is re- 
covered when decoherence is considered. This assertion was 
extensively tested numerically 0, E3, 13 l47ll mainly for the 
linear entropy (closely related to the purity) which is a lower 
bound of the von Neumann entropy. Other quantities, like 
the Loschmidt echo|5] which also display a noise indepen- 
dent Lyapunov decay, have also become of interest recently, 
especially in the context of quantum information processing 
and computing. Besides the linear entropy, in this section we 
study the asymptotic behavior of the autocorrelation function 
and the Loschmidt echo. 

For purely chaotic systems, after the initial spread gov- 
erned by the Lyapunov exponent, a state p evolved n times 
approaches asymptotically p^ = I/N and all time depen- 
dent quantities saturate to a constant value. The rate at which 
these quantities saturate is given by the largest eigenvalue, in 
modulus, smaller than one. Since, according to 11811 the spec- 
trum of L e approaches that of C e , then the universality of these 
decays can also be used to characterize quantum chaos. 

To display the decay towards p^ we subtract it from the 
initial state. Thus given an arbitrary state p, we define 

Po = P-jj, (42) 

were it is clear that Tr(po) = 0. Thus in all computations in- 
stead of evolving an initial state, we evolve an initial traceless 
pseudo-state such as the one defined in Eq. J42L orthogonal to 
Poo. Thus, we study how the distance between the initial state 
and the equilibrium state evolves. For example, for the lin- 
ear entropy, after the initial Lyapunov behavior, which ends 
at about the Ehrenfest time (tie ~ In N), instead of satura- 
tion to the equilibrium state p^, we expect to get an unbound 
growth which represents how this distance decreases exponen- 
tialy, and the exponent is proportional to |Ai|. 

Assuming for simplicity that all the eigenvalues are nonde- 
generate, and that Ri, i = 0, . . . , N 2 — 1, are the right eigen- 
functions (see App. then the expansion of po in terms of 
Ri is 

p = J2nRi, (43) 

were r» = Tr(Ljpo) and Li is the left eigenfunction. The 
pseudo state p evolved n times, is given by 

p n = L n e p = Y,n>%Ri. (44) 

If the eigenvalues are ordered decreasingly, according to 1 > 
I Ai| > | As j > • • • > Ajv 2 -i, then po is a sum of exponentially 
decaying modes. Suppose that Ai is real 3 , then it is clear from 
Eq. that 

p„ -> nX^Rt (45) 



3 In all the numerical simulations made, this was indeed the case. 
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FIG. 3: Quantum-classical correspondence for the noisy propagator. The top row shows repeated applications of the Perron-Frobenius operator 
of the perturbed Arnold cat of Eq. 1641 . to an initial classical (position) state. The bottom row shows the Husimi representation of po, . . . , p%, 
where po is a position eigenstate (N = 150, e = 0.2, k = 0.02). 



as n — * oo. Hence the asymptotic decay to the uniform den- 
sity is ruled by Ai. As a consequence any quantity which 
depends explicitly on p„ shows an exponential decay. Such is 
the case for the autocorrelation function 



C(n) = Tr(/3j/5 n ) 
From Eq. ( 145 \ we get, for large n, 

C( n )^\r 1 \ 2 X n 1 +. 



(46) 



(47) 



where we used the fact that Tr(R[Ri) = 1. If Ai is complex 
then 

and C(n) oscillates around A" (oscillation also appears if, for 
example, | A2 1 ~ |Ai|). Similarly, we can see that the linear 
entropy 



5„ = -ln[Tr(^)] 



(48) 



grows linearly with 2n. Once again, using Eq. i45l . the linear 
entropy for large n is 



In 



1 1 Ai ^(^fli) 



= -2nln[|Ai| 



constants. 



(49) 



Recently the Loschmidt echo has been extensively 
studied|5] especially in the context of fidelity decay in quan- 
tum algorithms 1 6j. The definition of the echo is 



M(t) = |(V>(0) |e^( H+s )*e'^ fft |^(0))| 2 



(50) 



which is the return probability of a state evolved forward a 
time t with a Hamiltonian H and backward with a slightly 
perturbed Hamiltonian H + S. It can also be viewed as the 
overlap between two states evolved forward with slightly dif- 
ferent Hamiltonians. Then M is just a measure of how fast the 
two states "separate". Most works focus on short times where 
several "universal" regimes have been identified. In particu- 
lar noise independent Lyapunov decay is observed for chaotic 
systems. 



In terms of the density operator, and discrete time systems, 
the Loschmidt echo after n steps is 



M n =Tr[p'„p n ]=Trf(U't) n (p )\J n (M 



(51) 



Where the prime represents a slight difference in the map. If 
the propagation occurs in a noisy environment, characterized 
by D e , it is natural to define the echo as 



A/„(e) =Tr 



(i.;t) n (p )L £ (po) 



(52) 



were Eq. ( 15 11 is recovered by making e = 0. 

Following the same arguments used for the autocorrelation 
function and for the linear entropy it can be shown that asymp- 
totically 



In [M n ]~n[ln(| Ai|) +ln(|Ai|)] 
Notice that Schwartz inequality implies that 



Tr [#,/>„] < \/Tr {fQ 



Tr 



(fir, 



(53) 



(54) 



Taking the natural logarithm of the expression above we get 



ln[M„] < I ( In ( Tr 



(Pr, 



In Tr 



(pn)' 



ln[M„] < --(S' n + S n ). 



(55) 



So we can see that the decay of the Loschmidt echo is bounded 
by the negative value of the average between the linear entropy 
of the original system and the perturbed one (see Fig. 4 in 1 8]). 

These three examples illustrate the fact that in the regime 
where the leading spectrum of L £ and C e coincide. We then 
expect all time dependent quantities to decay asymptotically 
with classical decay rates. 



B. Leading spectrum. Dynamics approach. 

In this section we describe the method used in 0] to com- 
pute the relevant eigenvalues of the coarse grained propagator. 
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This method works well for hyperbolic automorphisms of T 2 
because the nontrivial spectrum of the propagator lies entirely 
inside the unit circle for all values of e. The existence of a gap 
between 1 and Ai is crucial. 

In any complete basis, a superoperator such as L £ acting 
on Hn-2 has associated an iV 2 x iV 2 dimensional matrix. For 
small N this represents no setback. However, in order to es- 
tablish a relationship between quantum and classical we need 
to consider the semiclassical limit N — > oo and the diagonal- 
ization becomes unmanageable. To overcome this problem, 
we use an approach which takes advantage of the dynamics 
of the map to compute an approximation of the most relevant 
part of the spectrum by reducing sensibly the size of the eigen- 
value equation. 

Following l43l 14411 we construct two sets T, B € TLn 2 
which are explicitly adapted to the dynamics of the map 4 . 
Let po be an arbitrary initial density in Hn2, which for con- 
venience we choose it to be a pure state (projected onto some 
space). Then, by repeated application of L e we generate 



T = {fa, ft,..., ft, 

& = {P0,P{, ■ ■ ■ ,Pn, ■ 



where 



Pn = U(Pn-l) = L "(Pu) 

pi = Li(p° n _ 1 ) = o;(p ). 



(56) 
(57) 



(58) 
(59) 



Notice that L\ is the back-step propagator. Therefore, if 
the dynamics is chaotic, p" and p s n are increasingly smooth 
along the unstable and stable (classical) directions respec- 
tively. Thus they reflect the expected behavior of the left and 
right eigenfunctions of L e (see Fig.[3j. 

Using the bra-ket notation described in App. |A| we now 
construct the matrix 



'1,3 



(pt\Le\pJ ) = (Pt\ U e (pV) ) = (pj|$ ), (60) 



where (pf | = ( \J e (po)\ = (po\V t . Then we build the matrix 
of overlaps between elements of T and B, 

0« = (/a?|p?). (61) 
Notice that the structure of the matrices is very simple 

(p?|p?) = (poKlpf) = (po\K(pj)) 



(Po\p] +l ) 
( Pi+j I Pa ) 



(LI (p!)\po) 



(62) 



We remark that po G {poo}^- Because by construction L e is 
trace preserving, successive applications on an arbitrary po re- 
main in {poo} X and therefore the eigenvalue 1 is explicitly ex- 
cluded from our calculations. Moreover, the matrix elements 



4 See Florido, et al. [ 45] for a rigorous review on numerical methods that 
can be used to find RP resonances. The method used in l44ll . as well as its 
limitations, is analized there. 




FIG. 4: Plot of the matrix elements Oij = L £i(j ^ , where j+i — n. 
They are closely related to the autocorrelation C(n) = (po\Pn), 
Exponential decay is observed. The initial state is po = | 0, ) ( 0, | 
where | 0, ) is the coherent state centered at (0, 0), which is a fixed 
point of the map. 



in j60t and d6Udecav very rapidly, providing a natural cutoff 
n max to the sets T and B. In App. [5] we show that an approx- 
imation of the n max leading eigenvalues of L e arises from the 
solution of 



Det[[L e 



%3 



(63) 



i,j = 0,n n 



1, This method resembles the Lanczos itera- 



tion method 1 49] that uses Krylov matrices. 

The combination of small matrix computations plus a 
strong dependence on the dynamics makes this method a very 
efficient tool to get an approximation of the leading spectrum 
of L e for chaotic maps. 

Even when some of the main advantages of this method are 
evident (reduced size, leading spectrum and spectral decom- 
position, etc.), some drawbacks should be pointed out. When 
the classical system is nearly integrable some resonances can 
remain close to the unit circle and become unitary in the e — > 
limit and therefore convergence of the method with small ma- 
trices becomes problematic. Moreover, in that case there is 
a strong dependence in the initial state po. If it lies in a reg- 
ular island it will not explore all phase space. On the other 
hand, if initialized in the chaotic region it will only explore 
the chaotic sea, leaving out the regular tori. As a consequence 
some part of the relevant spectrum is inevitably lost. There- 
fore, the method is useful when the classical dynamics is fully 
chaotic. 



V. NUMERICAL RESULTS 

To illustrate our approach we utilize the Arnold cat map 
23l l with a small sinusoidal perturbation. The map is 



p' = p + q - 2nk sin[27rg] 

q'= q+p' + 27rfcsin[27rp'] 1 ' 



(64) 



where k is the small perturbation parameter. The map has a 
Lyapounov exponent which is almost independent of the value 
of k and equal to A = ln[(3 + v5)/2]. On the other hand the 
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FIG. 5: Leading spectrum of L e for different values of e and N. If Ai is the i-th eigenvalue, then log A; = log(ri) + i(f>i (where n = |A»|) 
and the coordinates in the plots are (</>, — log(r)). The ranges of the axes are cf> £ [— n, ir] and — log(r) £ [0, 6]. The map is the PAC with 
k — 0.02 and the matrix was truncated to dim= 12. 



Ruelle resonances (computed numerically) are very sensitive 
to it. Thus it is the ideal model to test the asymptotic results, 
independently of the short time Lyapunov regime. The map 
is a composition of two nonlinear shears and therefore it is 
easily quantized as a product of two noncommuting unitary 
operators. The explicit expression in the mixed representation 



The other advantage of using a map of this type is that the 
propagation both of pure states and of density matrices can 
be done by fast Fourier techniques, thus allowing relatively 
large Hilbert spaces with reasonable CPU times. The minor 
disadvantage is that the quantization for this particular map is 
only valid for even values of N 



[p\U\q] 



exp < i 



2vr 

TV 



qp 



P_ 

2 



cxp |27riVfc(cos[27n7/7V] + cos[27rp/7V]) j (65) 
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A. Spectrum 

In Ref. H| we have performed the classical calculation of 
resonances and shown that the quantum and classical leading 
spectra coincide. Here we take a slightly different approach 
and just compute the quantum spectra for a range of e and 
N values, as shown in Fig. [5] Observe that there is an ex- 
tended region where the spectrum is stable and independent 
of those parameters, signifying that the eigenvalues are prop- 
erties exclusive of the map, and therefore coinciding with the 
classical resonances. It is clear from this figure that the limits 
e — > and N — ► oo cannot be independent. In fact, at fixed N 
the limit e — > restores unitarity and the spectrum returns to 
the unit circle. Therefore, e must decrease as a certain func- 



// 


m 




/XT' 

i 




FIG. 6: Top (bottom) row shows the first 4 right (left) eigenfunctions 
showing the unstable (stable) manifolds for the quantum PAC with 
N = 100, e = 0.3, k = 0.02, and matrices truncated to dim= 12. 

tion of N. An optimal relationship between N and e is yet 
to be established but cannot be inferred from our limited data. 
However, in our range of values a dependence like e ~ 1/ ^/N 
seems suitable. 

The method described in Sect. IIV Bl also provides approxi- 
mations to the eigenfunctions of L e corresponding to the lead- 
ing eigenvalues. Inside the safe region (see Fig. [3 of N and e 
we were able to reconstruct at least 8 eigenfunctions success- 
fully with matrices of dimension of order 12. The accuracy of 
these eigenfunctions was checked by evaluating the orthogo- 
nality properties in Eq. 1381 and by computing the overlaps 



-{U\L t \R 3 ) 



(66) 



A plot of the absolute value of the Husimi representation for 
the first four eigenfunctions can be seen in Fig. [6] As was 
expected, the right (left) eigenfunctions corresponding to in- 
variant densities of the propagator is smooth along the classi- 
cal unstable (stable) manifold of the corresponding map. The 
right (left) eigenfunctions are not uniform along the unstable 
(stable) manifold showing pronounced peaks at the position of 
short periodic points. We intend to make a systematic analysis 
of this connection in a future work. 



Loschmidt echo for the PAC. In Fig. we see the growth of 
C n = — \n[C(n)] and the growth of S n for the perturbed 
Arnold cat defined in Eq. (I64> in Sect. [V] In both cases there 
are two well defined regimes. Initially both grow with the 
slope determined by the Lyapunov exponent of the map. For 
the PAC the Lyapunov is essentially the same for a wide range 
of perturbations. On the other hand, the Ruelle resonances de- 
pend strongly on the perturbation 5 . Taking as initial density a 
traceless pseudo-state (see Eq. (1421 ). time evolution of quanti- 
ties show how the state approaches uniformity exponentially, 
with a rate given by the largest RP resonance. We observe that, 
after the Lyapunov regime (around the Eherenfest time tie 6 ), 
the slope of the growth of C n is given by In | \\ | whereas the 
slope of S n is given by 2 In | Ai | as predicted. This factor two 
arises from the square in the definition of S n and is clearly 
seen in Fig.0 The solid lines represent these two slopes and 
were obtained by computing Ai using the method described 
in Sec UVBl 




FIG. 7: Purity decay and Correlation decay for the PAC (Eq. <64t ) 
with N = 450, e = 0.05, k = 0.005, initial pseudo-state p = 
| 0, 0)(0, — I/N, where j 0, ) is a coherent state centered at 
(0,0). The inset shows the evolution of S n for p — | 0, } ( 0, | and 
how it saturates to the constant value In N. 



In order to show the universality of the decay of the linear 
entropy and the Loschmidt echo, in terms of classical quanti- 
ties, in Fig.[8]we show S n and ln(M„(e)) for various values 
of the parameter e. The linear entropy is simply the negative 
logarithm of the purity Tr(p^). When e ~ the purity is con- 
served and equal to one, so the linear entropy does not grow. 
However when e ^ the purity will decayat a rate propor- 
tional to e. At one point, as predicted in [3], the growth of 
the linear entropy saturates and no longer depends on e. Since 
we evolved a traceless pseudo-state, with no component on 
the uniform density, after the Lyapunov growth the Ruelle- 
Pollicot regime appears. In the same way as for the entropy, 
for small values of e the asymptotic decay rate is e dependent 
but it saturates when rate determined by the first Ruelle res- 
onance is attained. This phenomenon can clearly be seen in 



B. Asymptotic decay 



In this section we study numerically the asymptotic behav- 
ior of the autocorrelation function, the linear entropy and the 



5 see Fig. 3 in @]. 

6 In Fig.QjV = 450 so n E ~ lniV = 6.11. 
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Fig- ID (left). In the right panel we display the echo and illus- 
trate exactly the same feature. 

VI. CONCLUSIONS 

We have developed a method to study numerically the spec- 
tral properties of open quantum maps on the torus. The 
method is particularly well suited to chaotic maps and pro- 
vides reliable eigenvalues and eigenfunctions. The noise 
model that we implemented utilizes phase space translations 
as Rraus operators and is equivalent to coarse graining quan- 
tum Markovian master equations. Therefore it brings out clas- 
sical properties of the map and we have shown that these prop- 
erties are reflected in the asymptotic decay of several quanti- 
ties. The same methods can be used to study other noise mod- 
els in the context of quantum information theory, if one thinks 
of the quantum map as an algorithm to be implemented and 
the noise as the error source present in any implementation. 
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APPENDIX A: ADJOINT AND LINEAR ACTION 

Let Hat be a complex Hilbert space of dimension N. The 
space of linear operators acting on Hn is called Liouville 
space £ = Elements in £ are usually represented 

by N x N dimensional complex matrices. However, given 



A, B £ H-n 2 then the "canonical" inner product, which in- 
duces the norm, is 

(A,B) =Tr(i f B). (Al) 

Thus H.N2 is a Hilbert space itself. Now, superoperators are 
a subset of the space of linear operators acting on H.n 2 ■ We 
introduce a bra-ket notation to simplify inner product expres- 
sions but also to distinguish the two types of decompositions 
we use for superopertors. Let A, B £ H-n 2 then the action of 
a superoperator S : TLn 2 —> "Hn 2 can t> e written as 

B = S(l) (A2) 

or as 

\B)=S\A) (A3) 

indistinctly. The adjoint, in the bra-ket form, is defined as 
usual by 

(A\S(B)) = (A\S\B) = (S\A)\B), (A4) 

which settles that (A\S = (S (A)\. Summarizing, 

A = \A) (A,-) = (A\ 

(A, B) = (A\B) S(i) = S\A) ' k ; 

One way to think about it (not absolutely necessary but help- 
ful) is to think of A as an operator, or matrix, in an operator 
space, acting on vectors, and | A ) as a vector in a vector space, 
acted on by superoperators. 

Now, a completely positive superoperator has a Rraus op- 
erator sum representation. Suppose S is a completely positive 
supeorperator then there exist a set of operators {M m }^ = q 1 S 
Hn2, such that 

S = ^M M ®Mt (A6) 
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are the Kraus operators. Without loss of generality, if the 
number of operators is smaller than A^ 2 we can always com- 
plete the set with zeros. The adjoint action of S on an operator 
A is defined through the Kraus representation suitable for the 
case of Eq. \A2l 



(A7) 



Eq. JA7> determines how the Kronecker product symbol <E> 
should be interpreted throughout this work. 

On the other hand, a superoperator S can be written as an 
expansion of spectral projectors. Let Ri and Lj be right and 
left eigenoperators of L respectively, such that 

S t (4) = A 4, z = l,...,iV 2 , (A8) 

S (£,) = AJL, 

and assume for simplicity that A, are nondegenerate. Then the 
spectral projectors are 

R i TT(Ll-) = \R i )(L i \, (A9) 
and the spectral decomposition is given by, 

s = J2\ R M L i\ < A1 °) 



(All) 



Therefore, given the spectral decomposition, the two ways of 
expressing the action of S on A are 

S| A ) = ]T I Ri ) A, ( Li\ A ) ee M^AMl = S(i) 

i i 

(A12) 

In more general terms Caves 1 39] identifies and describes the 
two different ways a superoperator acts as ordinary action 
(IA2> and left-right action (IA3I . This provides two distinct 
decompositions of the same superoperator. 



APPENDIX B: LEADING EIGENVALUES OF A LARGE 
MATRIX 

In this section we describe in a general way the method 
used to compute the leading eigenvalues of the superoperator 
L e in section I1VBI It is based on the Lanczos power itera- 
tion method|49] but was inspired by a recent work by Blum 
and Agam |44|. This method is useful when only a few of the 
largest (in modulus) eigenvalues is needed and also, since it 
deals with large matrices, when there is an efficient subrou- 
tine to implement the matrix-vector product but there is no 



need to store the whole matrix in an array variable. Moreover, 
convergence and accuracy depends strongly on the distance 
part of the spectrum one wants to calculate and the part to be 
neglected. 

In this work we don't address the question of the estimation 
of errors. 

Proposition 1 . Suppose A is a large, sparse matrix in C" x ™ 
and assume each of its eigenvalues \ has multiplicity one 
and that 1 > |Ao| > |Ai| . . . > |A n _i|. Suppose {/j}™_ and 
\riYl = Q are ^ ne corresponding left and right eigenvectors 



Ari = \ l r t 
A lj = Xjlj, 



(Bl) 
(B2) 



and let uq G C" be a vector such that 

|(Zj,uo)| > and \(n,u )\ > Vi < k (B3) 

for some k < n, where ( , ) represents as usual the inner 
product. Then the first k eigenvalues can be estimated from 
the reduced (k x k) eigenvalue equation 



Det 



K T (A\u ,k)A1C(A,uo,k) 

—zK r (A\u^, k)K(A, uq, k) 



= (B4) 



where K.(A, uq, k) is the Krylov matrix whose columns are 
the iterates of uq, 



K(A,u ,k) = u ,Au ,A 2 u ,...,A k u 
and T as usual denotes matrix transposition. 



(B5) 



Proof. We sketch a rather straightforward proof (though per- 
haps not entirely rigorous). The sets {Ii}™Zq and {ri}^T of 
left and right eigenvectors of A are complete and, they can be 
normalized according to 



(h,r 3 



(B6) 



Therefore there exist two distinct expansions of uo 

n-l 

u = '^a i r i (B7) 

i=0 
n-l 

(B8) 



i=0 
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In terms of these expansions we obtain 



K T (A\u ,k)K(A,u ,k) = 



EiA\* h 



(B9) 



which yields 
7C T (^4.^, mo, k)JC(A, uq, k) 

and similarly 

K T (A\u Q ,k)AK(A,u ,k) 



Thus, Eq. ( IB4t can be re-written as 



id 

i 

-- ^A^A^+^f .(Bll) 



If all the conditions of the proposition are met, this equation is 
equivalent to the original full eigenvalue equation. Now, since 
Xi are ordered by defeasing modulus and assuming that the 
eigenvalues accumulate around zero, leaving only a few, say k 
of them, with significant modulus (as is the case for the maps 
studied in section [^Hi} then we can neglect the contribution 
of the last n — k terms in the sum 7 . Thus Eq. JB 12i is just the 
determinant of the product of three k x k square matrices 



Det 

7 Although they can 



iq-copiputed, in this work we don't provide estimations 



0. 



(B12) 



Det[A f HA] = Det[A f ] Det[S] Det[A] = (B13) 



of the eiTors due to this truncation. 



where 



\1 A 



A 
Ai 

k-l 



Af 



v fc-i 



a^\ 



/ a p (X - z) 
' aiPx{Xi-z) 



\ 



(B14) 



a fc _i/3 fc _i(A fc _i - z) J 



The matrix A is a Vandermonde matrix. The determinant of a 
Vandermonde matrix Afc(Ao, . . . , Afc-i) is given by 

Det[A fc (A ,...,A fc _ 1 )] = J] (Xj-Xi). (B15) 

0<i<j<k-l 

From equation Eq. JB15i it can be readily seen that if the spec- 
trum of A is non-degenerate then A is invertible. Moreover, 
the structure of A determines k because in the limit of k "too 
large", A is singular, at least to within computing precision. 
So, using properties of the determinant in the secular equation 
Eq. dBT3l we get 

Det[S]=J]a^(A^-z)=0. (B16) 

Since, from the hypothesis, ^ then Eq. ( IB 16i yields 

the desired solution, i.e. the first k eigenvalues of A. □ 



In practice, the usefulness of the method depends upon the 
gap (1 — | Ai|), because it determines how fast the terms of the 
sum in Eq. ( IB 1 31 decay. 

In section HV Bl the span of the sets T and B are just the 
Krylov spaces l49ll of L c and \J e , and using the present notation 
Eq. (JS^J is 



Det 



1^ (Lj, PO, «max)L e /C(L e , pQ, n max ) 

-zX: t (L^po,"max)/C(L e , / 5 ,n max ) 



(B17) 



In analogy with Eq. dB 1 3 i . The efficiency of this method 
depends strongly on the spectrum configuration. The case 
of the coarse-grained propagator of hyperbolic maps on the 
torus 1 13 M 811 is particularly favorable because of the signifi- 
cant gap between 1 and Ai and because is an accumulation 
point, so a large number of resonances can be discarded and 
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the size of the matrices is reduced dramatically. 
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